ICM 2002 • Vol. Ill • 1-3 



Shock Waves 

Tai-Ping Liu* 



Abstract 

Shock wave theory was first studied for gas dynamics, for which shocks 
appear as compression waves. A shock wave is characterized as a sharp tran- 
sition, even discontinuity in the flow. In fact, shocks appear in many different 
physical situation and represent strong nonlinearity of the physical processes. 
Important progresses have been made on shock wave theory in recent years. 
We will survey the topics for which much more remain to be made. These 
include the effects of reactions, dissipations and relaxation, shock waves for in- 
teracting particles and Boltzmann equation, and multi-dimensional gas flows. 
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1. Introduction 

The most basic equations for shock wave theory are the systems of hyperbolic 
conservation laws 

«t + V a • /(«) = 0, 

where x 6 R m is the space variables and u € R n is the basic dependence variables. 
Such a system represents basic physical model for which u = u(x,t) is the density 
of conserved physical quantities and the flux f(u) is assumed to be a function of u. 
More complete system of partial differential equations takes the form 

u t + V x - f(u) = V x ■ (B(u, e)V x u) + g(u, x, t), 

with B(u,e) the viscosity matrix and e the viscosity parameters, and g(u,x,t) 
the sources. Other evolutionary equations which carry shock waves include the 
interacting particles system, Boltzmann equation, and discrete systems. Discrete 
systems can appear as difference approximations to hyperbolic conservation laws. 
In all these systems, shock waves yield rich phenomena and also present serious 
mathematical difficulties due to their strong nonlinear character. 
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2. Hyperbolic conservation laws 

Much has been done for hyperbolic conservation laws in one space dimension 

u t + f(u) x = 0, x G R 1 , 

see ^B], ID; HD, and the article by Bressan in this volume. Because the solutions 
in general contain discontinuous shock waves, the system provides impetus for the 
introduction of new ideas, such as the Glimm functional, and is a good testing 
ground for new techniques, such as the theory of compensated compactness, in 
nonlinear analysis. 

3. Viscous conservation laws 

Physical models of the form of viscous conservation laws are not uniformly 
parabolic, but hyperbolic-parabolic. Basic study of the dissipation of solutions for 
such a system has been done using the energy method, see [S] . Study of nonlinear 
waves for these systems has been initiated, [H], |14j . however, much more remains to 
be done. The difficulty lies in the nonlinear couplings due to both the nonlinearity 
of the flux f(u), which is the topic of consideration for hyperbolic conservation 
laws, as well as that of the viscosity matrix B(u, e). For instance, the study of zero 
dissipation limit e — > 0, see Bressan's article, has been done only for the artificial 
viscosity matrix B(u,e) = el. 

4. Conservation laws with sources 

Sources added to conservation laws may represents geometric effects, chemical 
reactions, or relaxation effects. Thus there should be no unified theory for it. When 
the source represents the geometric effects, such the multi-dimensional spherical 
waves , hyperbolic conservation laws takes the form 

^ m 

Ut + f(u) r = h(u), r 2 = y^(z») 2 . 

r i=l 

There has stablizing and destablizing effects, such as in the nozzle flows, 8 . The 
chemical effects occur in the combustions. There is complicated, still mostly not 
understood, phenomena on the rich behaviour of combustions. One interesting 
problem is the transition from the detonations to deflagrations, where the com- 
bined effects of dissipation, compression and chemical energy gives rise to new wave 
behaviour. Viscous effects are important on the qualitative behaviour of nonlinear 
waves when the hyperbolic system is not strictly hyperbolic, see Relaxation, 
such as for the kinetic models and thermal non-equilibrium in general, is interesting 
because of the rich coupling of dissipation, dispersion and hyperbolicity, 
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Conservative finite differences to the hyperbolic conservation laws, in one space 
dimension, take the form: 

- u n (x) = ^(F[u n ](x + ^) - F[u n ](x ^)). 

It has been shown for a class of two conservation laws and dissipative schemes, such 
as Lax-Friedrichs and Godunov scheme that the numerical solutions converge to the 
exact solutions of the conservation laws, [5]- On the other hand, qualitative studies 
on the nonlinear waves for difference schemes indicate rich behaviour. In particular, 
the shock waves depend sensitively on the its C-F-L speeds, |E|> |1^| - 

6. Multi-dimensional gas flows 

Shock wave theory originated from the study of the Euler equations in gas 
dynamics. The classical book [2] is still important and mostly updated on multi- 
dimensional gas flows. Because of its great difficulty, the study of multi-dimensional 
shocks concentrate on flows with certain self-similarity property. One such problem 
is the Riemann problem, with initial value consisting of finite many constant states. 
In that case, the solutions are function of x/t, not general function of (x,t), see 
0. See also |TB] for other self-similar solutions. However, unlike single space case, 
multi-dimensional Riemann solutions do not represent general scattering data, and 
are quite difficult to study. It is more feasible to consider flows with shocks and 
solid boundary, e.g. pQ 

7. Boltzmann equation 

The Boltzmann equation 

ft + &xf = Q(f,f) 

contains much more information than the gas dynamics equations. Nevertheless, 
the shock waves for all these equations have the same Rankine-Hugoniot relation at 
the far states. The difference is on the transition layer. There is beginning an effort 
to make use of the techniques for the conservation laws to study the Boltzmann 
shocks, [El- This is a line of research quite different from the intensive current 
efforts on the incompressible limits of the Boltzmann equation. 

8. Interacting particle systems 

Interacting particle systems is even closer to the first physical principles than 
the Boltzmann equation. There is the long-standing problem, the Zermelo paradox, 
in passing from the reversible particle systems to the irreversible systems such as the 
Boltzmann equation, and the Euler equations with shocks. This is fine for particle 
system with random noises. However, except for scalar models, so far the derivation 
of Euler equations from the particle systems has been done only for solutions with 
no shocks, |17|. 
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